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Abstract. We investigate the results of recently published experiments on the quantum behavior 
of Josephson circuits in terms of the classical modelling based on the resistively and capacitively- 
shunted (RCSJ) junction model. Our analysis shows evidence for a close analogy between the 
nonlinear behavior of a pulsed microwave-driven Josephson junction at low temperature and low 
dissipation and the experimental observations reported for the Josephson circuits. Specifically, we 
demonstrate that Rabi-oscillations, Ramsey-fringes, and spin-echo observations are not phenom- 
ena with a unique quantum interpretation. In fact, they are natural consequences of transients to 
phase-locking in classical nonlinear dynamics and can be observed in a purely classical model of 
a Josephson junction when the experimental recipe for the application of microwaves is followed 
and the experimental detection scheme followed. We therefore conclude that classical nonlinear 
dynamics can contribute to the understanding of relevant experimental observations of Josephson 
response to various microwave perturbations at very low temperature and low dissipation. 



1. Introduction 

In recent years the possibility of employing Josephson junctions circuitry as 
solid state elements for quantum information processing has been investigated by 
several groups and authors both at theoretical and experimental levels [1]. At tem- 
peratures low enough, operating the electronic Josephson devices in the regime 
where the thermal energy is smaller than the energy spacing between quantum me- 
chanical energy levels, interpretation of experimental results in terms of quantum 
response has been reported in a number of publications [2, 3, 4, 5, 6, 7, 8, 9, 10, 
11, 12]. Specific key observations that are being used to characterize the nature of 
a Josephson system at these extremely low temperatures are usually linked to the 
system response to a variety of pulsed microwave applications. The reason is at 
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least tri-fold. First, due to the nature of the Josephson effect, it is not easy to make 
direct observations of the expected quantum states in the anharmonic potential of 
a zero-mean voltage state Josephson system, which is the state of interest, so con- 
venient probes of the system are those where transitions from zero- voltage states to 
non-zero states can be induced with different probability depending on the nature 
of the initial zero- voltage state of the system. Second, the analogy between a single 
degree of freedom quantum Josephson system and a more traditional quantum sys- 
tem allows for useful characterization tools for Josephson applications of quantum 
mechanical features from, e.g., atomic systems, to infer the quantum behavior of 
the system from concepts such as Rabi oscillations [14], Ramsey fringes [15], and 
spin echo. Third, the control and manipulation of the state of the system is con- 
veniently directed through microwave frequencies that are commensurate with the 
anticipated energy level spacings in the zero-voltage Josephson state. The under- 
standing of the Josephson response to pulsed microwave applications is therefore 
critical for the interpretation of experimental observables. 

A classical nonlinear oscillator, which is subjected to pulsed time- varying per- 
turbations, may exhibit several phenomena that are relevant for microwave per- 
turbed Josephson systems. For example, an anharmonic oscillator provides an 
amplitude-dependent resonance frequency, which, in turn, relates resonance fre- 
quency directly to system energy [13, 16]. Also, the transient system behavior, 
due to an onset of resonant microwave application to a low-dissipation Joseph- 
son system, can give rise to low frequency amplitude and phase modulations in 
the system response [16, 17, 18, 19]. It has been shown that these phenomena in 
classical nonlinear dynamics are not unlike what one would expect from a compa- 
rable quantum system responding to a microwave signal, which is resonant with 
the intrinsic energy level spacing of the potential energy surface. In fact, the sim- 
plest possible classical Josephson system, a single Josephson junction, has been 
shown to exhibit phenomena of resonant escape from the anharmonic potential 
well [20, 21, 22], as well as transient slow modulations to phase-locking closely 
resembling experimental observations of Rabi-oscillations and Ramsey-fringes in 
Josephson circuits [16, 19, 23]. The correspondence is further substantiated by 
conducting numerical simulations of the classical Josephson model in accordance 
with the experimental procedures, and extracting the resulting data as it is done 
experimentally; namely as a statistical probability of escape from the Josephson 
potential well. Also, it has been observed that the classical model exhibits the phe- 
nomena of interest in the limits of low temperature and low dissipation (high-Q), 
consistent with reported experiments. With the demonstration of close classical 
analogies to Rabi-oscillations and Ramsey-fringes, which arc themselves closely 
related phenomena, we here extend the investigation of the classical model to also 
include a demonstration of the classical equivalent to the spin-echo observations 
[2, 4], since these represent yet another manifestation of the transient modula- 
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tions (Rabi-oscillations) in the approach to a phase-locked (excited) dynamical 
state. The work described here is primarily meant to demonstrate the phenom- 
ena in question in the most simple system possible; the single Josephson junction. 
However, while we will here limit ourselves to a descriptive treatment of a simple 
system, the phenomena are obvious and ubiquitous to a broad range of classical 
Josephson circuits, including one or more superconducting loops containing one 
or more Josephson junctions manipulated through magnetic fields imposed on the 
loops. 

2. The Model 

The normalized classical equation for a perturbed Josephson junction within 

(fi + Clip -\- sin (f = 1] + es{t) sin{uJst + 9s) + ep{t) -\- n{t) , 

where (p is the difference between the phases of the quantum mechanical wave 
functions defining the junction, rj represents the dc bias current, and Ss{t), cOg, and 
6s represent microwave amplitude, frequency, and phase, respectively. All currents 
are normalized to the critical Josephson current 7^., and time is measured in units of 
the inverse Josephson plasma frequency, wq"^, where luq'^ = 2eIc/hC = 2ttIc/^oC, 
C being the capacitance of the junction, and $o = h/2e = 2.07x10"^^ Wh is the 
flux quantum. Tunneling of quasiparticles is represented by the dissipative term, 
where a = fiujQ/2eRIc is given by the shunt resistance R, and the accompanying 
thermal fluctuations are deflned by the dissipation-fluctuation relationship [25] 

{n{t)) = (1) 
{n{h)n{t2)) = 2a-l—5{t2-h) = 2ae5{t2-h), (2) 

T being the thermodynamic temperature, Hj is the characteristic Josephson en- 
ergy Hj = Ich/2e, and O = ksT/Hj is thereby defined as the normalized temper- 
ature (or normalized thermal energy). In order to give orders of magnitudes for 
the typical realistic parameters we indicate that a good quality Josephson junction 
having a critical current of Jc = 5 and fabricated in the standard Nb-trilayer 
technology [26] may have R (usually identified with the subgap resistance) in the 
range (3 — Q)kVL , a capacitance C ~ bpF , vq = LOo/2Tr ~ 9GHz,a in the range 
(10-3 - 10"^), Hj ^ 16x10-22 J. The values of R an a in the specified interval 
depend on the quality of the junctions. 

A current pulse for probing the state of the system is represented by ep{t). Fig- 
ure 2 sketches the signaling details relevant for the observations of Rabi-oscillations 
(Fig. 2a), Ramsey-fringes (Fig. 2b), and spin-echo (Fig. 2c). The system is initially 
in its zero- voltage state (thermalized in the potential well illustrated in Figure lb), 
when a microwave is applied in a single interval or in a sequence of pulses, which 
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have the same reference phase Og ■ Due to the subsequent appHcation of the probe 
pulse Sp{t), the system variable (p may escape the potential well, and the ampli- 
tudes of ep{t) is parameterized such that probe pulse induced escape is likely if 
the system energy, relative to the potential well depth, is significant at the time 
immediately prior to the onset of ep{t), and unlikely otherwise. By conducting a 
large number of such simulations, each with identical parameters, except for Og, 
which is chosen from a uniform distribution in the interval ]0, 27r], and the ther- 
mal noise current, a statistical signature of the system energy due to the signaling 
strategy given by £sit) can be inferred from the resulting switching probability. 
This information can be directly obtained in simulations by calculating the defined 
normalized energy 

1 2 

H = -^ip -I- 1 - cos V9 - ry(/9 , 

where the minimum of the potential well can be represented by the energy 

Hq = 1 — a/I — 77^ — r] sm~^ rj . 

The above statistical procedure of monitoring the probability of escape is applied 
experimentally, since direct measurements of the energy content in a zero-average 
voltage state is not possible. Wc therefore mimic the experimental procedure in 
our simulations in order to make direct linkage between modeling and experimen- 
tal observations of switching probabilities. We finally notice two essential system 
parameters that are related to the choice of bias point 77; namely microwave fre- 
quency LOg and temperature Q. The microwave frequency should be chosen such 
that resonant states can be excited in the vicinity of iOg « cof = sjl — rf, and 
the thermal energy should be significantly smaller than the potential well depth 

sin ^rj — Trrj, ensuring that recorded switching events 
are not simply a reflection of thermal activation. 

3. Simulation Details and Results 

Following the procedure of reported experiments, we record the switching from 
the zero- voltage state ((0) = 0) as a result of applying the three different recipes 
sketched in Figure 2. We will, unless otherwise noted, throughout this presenta- 
tion use the following parameters: r] = 0.904706, which is equivalent to a linear 
resonance frequency cj/ = ^1 — r/^ = 0.652714, a = 10~^, which provides a high-Q 
resonator, and, when temperature is applied, = 2 x 10"*^ -C AU = 5.57 x 10"^. 

3.1. Rabi-Oscillations 

The first signaling procedure is shown in Figure 2a, where the Josephson junc- 
tion (Equation (1)) is initially resting in the potential well, when the microwave 
pulse is initiated. We use a microwave frequency LOg = ui, but we note that this 
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is not a requirement for the observations of tliis work [16]. Figure 3 shows the 
evolution for = of the junction phase ip{t) and the system energy H{t) — Hq 
for a microwave pulse with amphtude = 0.00217. Several important features 
can be extracted from this figure. First, we observe how the junction is driven 
into a phase-locked state by the microwave application. Second, we observe how 
the junction response is modulated by a low frequency envelope function that 
causes the system energy to oscillate around the steady-state phase-locked energy. 
Third, comparing two different durations of microwave pulses, we illustrate how 
the microwave duration inflTicnces the switching probability when the probe pulse 
is applied after the termination of the microwave field. Figures 3a and 3b show 
how the probe pulse catches the system in a low energy state, resulting in a "no- 
switch" count, and Figures 3c and 3d show the probe pulse catching the system 
in a high energy state, resulting in a "switch" count. Thus, by varying the du- 
ration of the microwave pTilsc, we can observe the slow modulation frequency of 
the transient to the phase-locked state in the zero- voltage state of the junction by 
observing the switching statistics by the probe pulse application. Such statistics is 
shown in Figure 4a for 6 = 2 x 10^"^, where the switching probability as a function 
of microwave duration is displayed for averages of 2500 possible switching events. 
Clearly, the probability exhibits the slow modulation seen in Figure 3, and these 
are the classical analogues to Rabi-oscillations [16, 19]. These oscillations have 
a reasonably well-defined frequency. The amplitude decays for several reasons, 
including the dissipation a, which sets a time scale for the transients to phase- 
locking [16], thermal effects, which provides decoherence to the evolution, and the 
random phase 9s, which dephases the different trajectories that are averaged to 
give Figure 4a. Based on Figure 4a, we can now extend the analogy to quantum 
mechanics by noticing that a microwave pulse of one Rabi-period tji = 2TT/i}fi, 
where ^Ir is the modulation (Rabi) frequency, leaves the system energy in roughly 
the same state as before the microwave pulse was applied. We therefore denote 
such a pulse a 27r-pulse. Similarly, we define a 7r-pulse to be one of duration tr/2, 
and a 7r/2-pulse represents a duration r/j/4, which is the duration it takes from 
the microwave onset to elevate the system to its steady state energy. 

3.2. Ramsey-Fringes 

Following the recipe outlined in, e.g., Ref. [2, 4] we can now apply a sequence of 
controlled microwave pulses, such as iUustrated in Figure 2b, where two ir/2 pulses 
arc applied with a delay before tlic probe pulse determines if the resulting state 
has enough energy to be perturbed into escape. Figure 5 illustrates, in analogy 
with Figure 3 for = 0, the two characteristic outcomes of temporal evolutions 
of the system dynamics for two different time intervals between the 7r/2-pulses. 
The first pulse elevates the system energy to the intermediate level, where it is 
left to evolve freely after the first pulse is terminated. During the time between 
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the two pulses, the microwave field and the junction oscillate at almost the same 
frequency, but with a small detuning of the mutual phase. Thus, this ballistic in- 
terval can produce very different outcomes when the second 7r/2-pulse is applied, 
depending on the time interval. If the time interval is such that the second pulse 
is initiated with a mutual phase difference to the junction similar to when the 
first pulse terminated, then the second pulse will continue to pump energy into 
the junction, elevating its energy to the maximum at the conclusion of the second 
pulse. This scenario is shown in Figures 5c and 5d, where the probe pulse triggers 
an escape event. However, if the second 7r/2-pulse is initiated such the the relative 
phase is opposite to what it was when the first pulse was terminated, then the 
second pulse will attenuate the energy of the junction, leaving it with very low en- 
ergy at the conclusion of the second pulse. This scenario is illustrated in Figures 
5a and 5b, where the probe pulse cannot induce an escape event. The switching 
probability as a function of the ballistic time interval between the two 7r/2-pulses 
is shown for = 2x10^'* and microwave frequency Ug = u)i in Figure 4b. The 
clearly visible oscillations in the switching probability are the classical equivalent 
of quantum mechanical Ramsey- fringes. As was the case for the Rabi-oscillations, 
the visible Ramsey- fringes have a reasonably well defined frequency Vl p and an at- 
tenuating amplitude, where the attenuation is due to a combination of dissipation, 
thermalization, and the averaging over trajectories of different random microwave 
phases Og- The close relationship between the classical observation shown here and 
quantum mechanical Ramsey-fringes can be further underlined by evaluating the 
Ramsey- fringe frequency f^i? as a function of frequency detuning form the har- 
monic resonance frequency oj;. This is shown in Figure 6, where black markers 
and associated error bars represent the measured Ramsey-fringe frequency. We 
clearly observe the characteristic unity-slope ("V" -shape indicated by the dashed 
lines) dependency on the frequency detuning from linear resonance [10]. An obvi- 
ous deviation from a naive expectation is that the microwave frequency for which 
Of = is not exactly the linear resonance. This is due to the anharmonicity of 
the potential as well as to the dissipation in the system. The Ramsey-fringe fre- 
quency fii;' is drifting slightly toward higher values as the 7r/2-pulse separation is 
increased. This is because higher oscillation amplitudes result in lower resonance 
frequencies. Thus, defining the Ramsey-fringe frequency from the first couple of 
oscillations (the ones with the best resolution) will make Q.f susceptible to the 
anharmonicity of the potential. This, in turn, will make the ballistic evolution 
between the two 7r/2-pulses commensurate with a slightly detuned microwave fre- 
quency such that fip = for a;^ < w^; which is what can be observed in Figure 6a. 
Figure 6b illustrates how increasing the dissipation parameter to a = 10^^ will 
provide a faster linearization of the resonance frequency of the junction, resulting 
in the expected shift in the characteristic point where $7^ = for ~ t*^;- 
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3.3. Spin-Echo 

We finally study the result of applying the signaling protocol shown in Figure 
2c. This protocol is similar to the one resulting in Ramscy-fringcs, but a vr-pulse 
is inserted between the two 7r/2-pulses in order to mimic a 180° phase shift in 
the modulation dynamics. In analogy with experimental observations of spin- 
echo in Josephson systems [2, 4] we choose a temporal separation between the 
two 7r/2-pulses, and observe the switching characteristics of the junction, when 
the probe pulse is applied, as a function of the temporal position of the 7r-pulse, 
which serves as a 180° phase-shift from which an "echo" can be detected. Since 
the system either switches or remains in the potential well after the application of 
the probe, wc can again demonstrate the typical dynamics by the two cases shown 
for = in Figure 7. Here we have chosen a total normalized time interval for 
the 7r/2-pulses of 2500 time units. The first case, which is displayed in Figures 7a 
and 7b, shows that the application of the 7r-pulse only has a marginal effect on 
the evolution of the phase and its energy. Since the second 7r/2-pulse is initiated 
out of phase with the junction, the subsequent probe pulse perturbation cannot 
make the system escape the potential well. In contrast. Figures 7c and 7d show 
that a slightly different timing in the application of the 7r-pulse will provide a 
resonant increase in the system oscillation amplitude and energy. Further, the 
phase-relationship between the microwave field and the junction is now such that 
the second 7r/2-pulse enhances the energy of the system. This results in an escape 
event when the probe pulse is applied. The probability of switching as a function 
of the timing of vr-pulse initiation is shown in Figure 4c for = 2 x 10~^ and 
ujs = 0.965cJi. We find very clear signatures of the oscillations in the switching 
probability, which can be interpreted as the phase-dependency of the "echo" that 
either leaves the system in high or low energy states. Thus, this is the classical 
Josephson analogue to quantum mechanical spin-echo. We again notice that the 
oscillations provide a certain frequency, which we have displayed in Figure 6a as 
open markers as a function of microwave frequency detuning. It is obvious that the 
frequency is closely related to the Ramsey-fringe frequency ftp, except near the 
point ^Ip = 0. This deviation can be directly associated with the drift in frequency 
due to anharmonicity that we mentioned above for the Ramsey-fringes, and the 
exact values of the spin-echo oscillations seem to depend strongly on which part of 
the curve (such as the one shown in Figure 4c) is used to calculate the frequency 
when the oscillation frequency is small. 



4. Conclusions 

Just like spectroscopic resonant multi-peak switching distributions can be ob- 
served and explained within a classical Josephson framework [20], so can the phe- 
nomena of Rabi-oscillations, Ramsey-fringes, and spin-echo be observed in the 
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purely classical Josephson model by following the recipe of the experimental ob- 
servations. We have here tried to outline this observation in the simplest possible 
system; namely the single Josephson junction. This is not a direct representation of 
specific experimental configurations, which are much more complex, although the 
parameter values chosen here have been inspired by several relevant experiments. 
Considering the simplicity of the single junction and the modeling approach in 
this presentation, the overall agreement between the classical results and experi- 
mental observations is remarkable. Not only do we observe the abovementioned 
phenomena in the simple classical model, but we also observe the magnitudes 
of, and relationships between, the various oscillations in the respective switching 
distributions to be quite reasonably aligned with experimental observations and 
expectations. We therefore submit that much in the interpretation of the ex- 
periments on microwave manipulations of Josephson systems at low temperature 
and low dissipation can be contributed from the well known and well documented 
classical Josephson dynamics. 
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Fig. 1: Sketch of the RCSJ model where I/Ic = f] + es{t) sin^ajst + 6) + £p{t) and 
N = n{t)Ic- Also shown is the Josepshon potential for rj = 0.904706. 
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Fig. 2: Signalling diagrams for simulating (a) Rabi-type oscillations (b) Ramsey- 
type fringes and (c) spin-echo-type oscillations. 
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Fig. 3: Direct simulation of the Rabi-type oscillation switching response. Joseph- 
son junction response to the sequential application of two 7r/2 microwave pulses 
followed by a probe field. Panels (a,b) show phase-difference and energy for a non- 
switching sequence. Panels (c,d) indicate a switching event. Switching and non- 
switching events trigger according to a randomly determined phase. Parameters 
were a = IQ'^, r] = 0.904706, = 2.17 x 10"=^, Us = ui = ^1 - r/2 = 0.652714, 
£p = 8.2 X 10-2, and 6 = 0. 



Classical analogs for Rabi- oscillations, Ramsey-fringes, and spin-echo in Josephson junctions 13 




Fig. 4: Switching probabilities for (a) Rabi-type oscillations, (b) Ramsey-type 
fringe and (c) Spin-echo-type simulations. Each point represents statistics of ~ 
2, 500 events at 9 = 2.00 x 10"''. The horizontal axis in panel (a) represents the 
duration At^ of the microwave pulse. In panel (b) the horizontal axis indicates 
the 7r/2-pulse separation (i4 — ti). Panel (c) has the period between the end 
of the first 7r/2-pulse and the beginning of the 7r-pulse - i.e., the 7r-pulse offset 
{t2—ti). In panels (a) & (b) the driving frequency is the linear resonance frequency 
oog = ijJi = yl — ry^ = 0.652714. In (c) the driving frequency is oog = 0.965a;;. The 
magnitude of the probe pulse in (a) & (b) Ep = 8.2 x 10~^. The probe pulse in 
panel (c) is £p = 7.5 x 10~^. Panel (c) reflects data from spin-echo-type simulations 
for t^ — ti = 2750. Common to all panels are the characteristic damping a = 10~^ 
and the microwave signal amplitude varepsilous = 2.17 x 10~^. Lines are added 
to aid the eye. 
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Fig. 5: Direct simulation of the Ramsey-type fringe switching response. Joseph- 
son junction response to the sequential application of two 7r/2 microwave pulses 
followed by a probe field. Panels (a,b) show phase-difference and energy for a non- 
switching sequence in the case of t^^ — ti = 400. Panels (c,d) indicate a switching 
event for -ti = 7 0. Parameters were a = lO""^, rj = 0.904706, Eg = 2.17 x 10"^, 
uj, = uji= ^l-r/2 = 0.652714, Sp = 8.2 x 10"^ and 9 = 0. 
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Fig. 6: Frequency response as a function of applied microwave frequency for two 
different dissipation parameters: (a) a = 10~^, and (b) a = 10""^. Filled dia- 
monds depict Ramscy-typc fringe frequencies, Jl^ , with parameters corresponding 
to Fig. 4. The open circles in (a) represent spin-echo-type frequencies with param- 
eters as in Fig. 5. 
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Fig. 7: Direct simulation of the spin-echo- type oscillation switching response. 
Josephson response to the application of two 7r/2 microwave pulses, with an in- 
tervening vr-pulsc, followed by a probe field. The delay between the 7r/2 pulses 
is t/^ — t\ = 2182. The vr-pulsc offset is for (a,b) (non-switching) t2 — t\ = 1000. 
For (c,d) (switching) t2 — t\ = 900. Parameters were a = 10-4, ^ = 0.904706, 
Es = 2.17 X 10-=^, £p = 7.5 X 10-2, and 9 = 0. 



